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Linear Polarization as a Probe of Gamma Ray Flaring Blazar Jets 
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We describe and present initial results from a Fermi cycle 2 program designed to monitor the behavior of the 
centimeter-band linear polarization and total flux density emitted by gamma-ray-bright blazars during flaring. 
The goal of the program is to identify changes in the magnetic field structure in the radio jet associated with 
gamma-ray flaring and ultimately to test whether gamma-ray flaring is associated with the onset of shocks in 
the radio jet. Light curves illustrating radio band variability patterns are shown for sample program sources. 
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1 . Overview of Project 



As part of a Fermi cycle 2 program, we are obtain- 
ing monitoring observations at three centimeter band 
frequencies (14.5, 8.0 and 4.8 GHz) with the Uni- 
versity of Michigan 26-m radio telescope (UMRAO) 
of the source-integrated total flux density and linear 
polarization of approximately 30 radio- and 7-ray- 
bright sources currently or expected to be in 7-ray- 
flare phase. The goal of our work is to identify jet 
conditions responsible for the generation of the 7-ray 
emission and to test the hypothesis that shocks play 
a role in the production of this emission as suggested 
by earlier work 0]. A causal connection between ac- 
tivity in the radio jet and T-ray flaring first proposed 
based on EGRET data Sj is supported by early 
Fermi results. For example, statistical studies of the 
sources in the 15 GHz MO J AVE VLB A sample 0| 
have found that the members detected during the first 
three months of Fermi operation 0] are more core- 
dominated, have higher brightness temperatures and 
Doppler factors, and are in a more active radio state 
than non-detected MO J AVE sources 0,0]. 

While there has been increasing evidence suggesting 
that the 7-ray emission arises near the radio core 0], 
a region believed to correspond to either the t—1 sur- 
face or to a standing shock 0] and not at a site near to 
the black hole/ accretion disk, a number of important 
questions remain. These include the trigger for the 
emission of 7-rays, the position of the emission site 
(upstream or downstream of the radio core) and the 
nature of the emitting region itself (standing shock, 
propagating shock, or 'blob' with a chaotic magnetic 
field where turbulence produces the very high energy 
electrons). A plausible scenario for the broadband 
features seen as outbursts or flares is that instabilities 
develop naturally within the collimated relativistic jet 
outflows giving rise to shocks [TO] - With the passage 
of a shock, there is a compression of the magnetic field 
within the emitting region and an increase in the de- 
gree of order of the magnetic field. The observational 
signature of such an event is a swing in the electric 
vector position angle (hereafter EVPA; an orientation 
orthogonal to the magnetic field direction in a trans- 
parent source) and an increase in the fractional linear 



Table I Program Sources 



0048-097 


NRAO 190 


0917+449 


1510-089 


2022-077 


0109+224 


0454-234 


1127-145 


1633+382 


BL Lac 


DA 55 


0528+134 


3C 273 


3C 345 


2201+171 


0215+015 


0716+714 


3C 279 


1717+178 


CTA 102 


3C 66A 


0727-115 


1308+326 


OT 081 


3C 454.3 


0235+164 


0805-077 


1502+106 


1849+670 




3C 84 


OJ 287 


1508-055 


1908-211 





polarization. 



2. The Source Sample 

Our program sources are listed in Table 1. In addi- 
tion to blazars, our source sample includes the radio 
galaxy 3C 84 (NGC 1275); 7-ray emission has been de- 
tected in this source which appears to be associated 
with compact regions within the radio jet [HI]. Selec- 
tion was based on the following criteria: 1) inclusion 
in the high confidence list in the LAT catalogue based 
on the first three months of Fermi's operation 0] , or 
identified as 7-ray flaring in subsequent Astronomers' 
Telegrams; 2) 15 GHz total flux density of at least 1 
Jy in early 2009- and 3) membership in the MOJAVE 
15 GHz survey [5j. While AGN light curves regularly 
show peaks and troughs that can last for months or 
longer, we expect these objects to be highly likely to 
exhibit bright flaring states during cycle 2. A mini- 
mum radio band total flux density of 1 Jy was chosen 
to ensure good signal-to-noise in the multifrequency 
polarimetry data. Typical fractional linear polariza- 
tions in blazars are of order a few percent, and the 
radio band spectra of our targets are characteristi- 
cally flat or inverted (14.5 GHz total flux density > 
4.8 GHz total flux density). The MOJAVE 15 GHz 
imaging data, typically obtained every few months or 
longer, will complement our single dish measurements 
which are more frequent but which integrate over the 
entire source area. About half of our sources are in the 
Boston University 43 GHz VLBA program, providing 
information on morphology and structural changes in 
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Figure 1: From bottom to top, daily averages of total flux density S, polarized flux P, and electric vector position angle 
X illustrating recent variability at radio band in the high redshift (z=1.839) QSO 1502+106 (OR f03). Observations at 
14.5, 8.0, and 4.8 GHz are denoted by crosses, circles, and triangles respectively. 



the inner jet. 

As new 7-ray-flaring radio-bright AGN are iden- 
tified, they are added to the UMRAO program if 
they meet our criteria for inclusion in the monitor- 
ing program; examples of newly added sources are 
NRAO 190, a bright EGRET source Q3 in which 
recent -Fermi-detected flaring was reported [l3j |. and 
3C 345 a well-known bright blazar surprisingly not 
detected by EGRET. Additionally, we are observing 
on a time-available basis a few sources which were ex- 
ceptionally bright during the EGRET era but which 
are currently relatively quiescent in both the radio and 
7-ray bands. Included in this group is NRAO 530 [14] . 
Several of the sources listed in Table 1 are new to our 
core program; others have been observed by UMRAO 
for many years. 



of 8.0 and 4.8 GHz for flaring sources. This sampling 
rate is matched to the expected variability timescale 
and duty cycle in the radio band based on historical 
measurements obtained over decades by UMRAO [l5[ ■ 
However, the sampling rate will be increased if war- 
ranted, e.g. during some phases in the BL Lac object 
0716+714 which has historically exhibited large am- 
plitude, very rapid radio band variations. Each daily 
measurement consists of a series of on-off measure- 
ments over a 30-40 minute time span; these measure- 
ments are interleaved with observations of a calibra- 
tor every 1-2 hours to monitor the antenna gain and 
pointing, and to verify the instrumental polarization. 
As a result of these requirements, only 20-25 program 
sources can be observed in each 24 hour run, and the 
most active will be selected from the sample. 



2.1. Observing Strategy 

Before selecting our sources, we carried out ex- 
ploratory radio band measurements of a larger group 
in 2009 March to evaluate the current radio band vari- 
ability state and to measure the flux levels for all po- 
tential program sources. Monitoring of our sample 
commenced in 2009 August. Throughout cycle 2 we 
expect to obtain two observations weekly at 14.5 GHz 
(the frequency at which the variability amplitude is 
highest in AGN), and one observation weekly at each 



3. Early Program Results 

The first stage of our program, to obtain data ex- 
hibiting resolved linear polarization variability tem- 
porally associated with time periods of 7-ray flaring, 
is in progress. During the first two months of the 
program we have detected several polarization events, 
and we are in the process of examining these and our 
subsequent data to identify characteristic patterns in 
the light curves. Resolved flares potentially suitable 
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Figure 2: From bottom to top, daily averages of total flux density, fractional linear polarization, and electric vector 
position angle at three frequencies illustrating the recent variability at radio band in 3C 279. Note the changes in 
polarization in 2009 April-May which are resolved in these single dish observations. The variations in this source 
during the mid 1980s were successfully modeled with a transverse shock [l7| . 



for detailed modeling have been identified in 0727-115, 
0805-077, 3C 279, 1502+106, and 1510-089. 

We are finding that the timescales of the outbursts 
in polarized flux and periods of ordered temporal 
changes in EVPA are relatively short, typically sev- 
eral weeks to a few months in duration. As an exam- 
ple, we show in Figure[T]the light curves for 1502+106 
(OR 103), a relatively new addition to our core pro- 
gram, which is highly active at 7-ray band. The MO- 
JAVE 15 GHz source structure in this source is rela- 
tively simple, and much of the polarized emission is 
associated with a single source region. A study of 7- 
ray flaring in August 2008 used two MOJAVE images 
separated by several months to identify a change in 
magnetic field orientation [l6l |. but the sampling of 
the imaging data did not permit tracing the change in 
structure. During the subsequent period covered by 
our monitoring data, there is a resolved outburst in 
polarized flux (middle panel) and a systematic swing 
in EVPA (top panel) which commenced in 2009 Febru- 
ary. These variations track at 14.5 and 8.0 GHz; but 
the data follow a different path at 4.8 GHz. The multi- 
frequency total flux density data (bottom panel) show 
a spectrum characteristic of a self-absorbed source; 
thus the emission at 4.8 GHz most likely arises in a 
somewhat different physical region of the source, fur- 
ther out from the inner region probed by the higher 



frequencies. In this flare the position angle swing and 
the development of the linearly polarized outburst are 
characteristic of a shock event. 

Figure [5] shows results for a structurally complex 
QSO, 3C 279, where the contributions from individ- 
ual core and jet components are blended in the total 
flux density light curve. During a series of events in 
the mid-1980s the variations apparent in the UMRAO 
data for this source were successfully modeled assum- 
ing a transverse shock [l7|; and the jet parameters 
derived in that work will be used as a starting point 
in the new analysis. 

Some sources exhibiting 7-ray flaring have not 
shown large amplitude changes in either total flux den- 
sity or in polarized flux in our observing band (e.g. 
NRAO 190). We do not know yet whether this is 
due to frequency-dependent time delays across bands, 
the masking of variability in the UMRAO source- 
integrated measurements due to competing, indepen- 
dently evolving source components, or whether more 
than one scenario is required to explain the origin of 
the high energy emission. 

Two very bright EGRET sources on the Fermi LAT 
list, 0528+134 and NRAO 530, have not exhibited 
high-amplitude flaring in either the 7-ray or in the ra- 
dio band since the launch of Fermi. We show the long 
term light curves for NRAO 530 in Figure [3] which il- 
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Figure 3: From bottom to top, long term monthly averages of the total flux density, polarized flux, and electric vector 
position angle in the QSO NRAO 530, a source which was exceptionally bright in both the 7-ray and in the radio band 
during the EGRET era in the 1990s, but which is currently exhibiting only low level variability in both of these bands. 



lustrates the relatively low flux level in the radio band 
since the launch of Fermi. The presence of low fluxes 
in both bands is consistent with the expected behav- 
ior assuming that the activity is broad band and that 
the same particles are responsible for the high and the 
low energy emission. 

In the next phase of our work, we will carry out 
shock modeling of resolved radio band events exhibit- 
ing the shock signature. This will allow us to set limits 
on the physical conditions in the radio jet during 7- 
ray flaring. A set of resolved radio band flares will 
be selected for detailed shock-in-jet modeling follow- 
ing the procedures in our previous analyses [l^, 53] 
which use multifrequency linear polarization and to- 
tal flux density observations as constraints. While our 
earlier work assumed a specific shock geometry (trans- 
verse shocks), the new models will employ a transfer 
code that admits arbitrary shock orientation, result- 
ing in extra degrees of freedom which must be con- 
strained. The complementary VLBA imaging data 
from the MOJAVE and BU programs will be used to 
limit the allowable range of orientations. The model- 
ing is expected to yield information on shock strength, 
the character of the jet flow, and the low energy cut- 
off in the spectrum of the radiating particles, all of 
which are key input parameters in jet emission mod- 
els designed to explain the origin of the 7-ray emission 
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